Dilepton-tagged jets in relativistic nucleus-nucleus collisions: A case study 
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We study the A + B — > £ + £~+ jet +X process in nucleus-nucleus collisions at relativistic energies. 
The dilepton as well as the jet will pass through the matter produced in such collisions. The recoiling 
dilepton will carry information about the kinematical features of the jet, and will thus prove to be 
a very effective tool in isolating in-medium effects such as energy-loss and fragmentation function 
modifications. We estimate the contributions due to correlated charm and bottom decay and we 
identify a window where they are small as compared to pairs from the NLO Drell-Yan process. 



I. INTRODUCTION 

The study of dilepton production through the Drcll- 
Yan process p], hi + hi — > £ + i~ + X has remained 
a useful tool to sound out concepts about the parton 
model, QCD, and the structure functions of hadrons. If 
the initial state partons have no transverse momentum, 



the lowest-order process qq — > 7* 



I produces a lep- 



ton pair with a net Qt=0- Experiments, however, show 
that the net transverse momentum of the dileptons pro- 
duced by the Drell-Yan process are of the order 1 GeV for 
a dilepton having a mass M ~ 10 GeV 2j. This fact out- 
lines the importance of the NLO Drell-Yan contributions. 
Also, one could assign an intrinsic spread to the trans- 
verse momenta of the partons. This momentum can have 
many origins, and a confinement argument along those 
lines only accounts for < q\ >~ (0.3 GeV) 2 Q- 

It is therefore clear that the dileptons acquire addi- 
tional transverse momenta from production mechanisms 
beyond leading order in perturbation theory: 



qg -> 97 



and qq — > 57* 



(1) 



These Compton and annihilation processes arc analogous 
to the ones responsible for the production of real photons 
where the recoiling final state quark or gluon balances 
the qr of the dilepton. An important point £j is that 
the dileptons thus produced are always accompanied by 
a recoiling quark or gluon. If the energy of the quark or 
the gluon is several tens of GeV, it will lead to a jet of 
hadrons in a narrow cone around the leading hadron in 
the jet. 

Now consider a collision of two heavy nuclei at rela- 
tivistic energies, which could lead to the formation of a 
quark gluon plasma (QGP). The dilepton and the jet pro- 
duced in the above process will pass through the plasma. 
The dileptons will not interact, while the jet will lose en- 
ergy through collisions and radiation in the plasma and 
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provide valuable information about these mechanisms in 
the medium. The dilepton can thus be used to tag the 
jet in the same manner that photon-tagged jet measure- 
ments have been suggested as a precise probe for the 
study of jet-quenching [jj. 

We recall that jet-quenching can manifest itself in vari- 
ous ways; there would be a suppression of hadrons having 
large transverse momenta, when compared to the results 
for pp collisions, and the fragmentation function in AA 
collision may be different from those seen in e + e~ or pp 
collisions, etc A precise determination of these ef- 
fects will require knowledge of the energy of the parton 
produced in the hard collision which fragmented into the 
jet. Note that the inclusive transverse momentum distri- 
bution of the hadrons will not permit us to deduce the 
value for the parton dE/dx unambiguously, as a hadron 
with a given p^, adron can arise from the fragmentation of 
any parton having p^ rton = p^? dlon /z, where p^ rton is 
the momentum of the parton at the time of fragmenta- 
tion, and z (< 1), is the fragmentation variable. Recall 
also that the vast amount of theoretical activity [(| in 
the last decade has yielded varying predictions for the 
dependence of dE/dx on the energy of the parton and 
the properties of the medium. A precise knowledge of 
the energy of the progenitor of the jet of hadrons will 
provide the information necessary to settle these issues. 

The very large background of hadrons in AA collisions 
makes it very difficult to characterize the jet through 
canonical jet algorithms which typically prescribe energy 
measurements in some cone around the jet-axis. To min- 
imize contributions from the large non-jet background, 
particles emitted beyond the cone or having an energy be- 
low some threshold are excluded in the jet measurement. 
This excluded energy, as well as the non-jet background 
energy, must be estimated and taken into account. Be- 
cause of the difficulty of making these corrections, photon 
tagging of the jet is a valuable alternative or complement 
to the jet energy measurement, as the transverse momen- 
tum of the photon (real or virtual) is equal to that of the 
high energy parton produced in the collision. 

In this work, a measurement of dilepton-tagged jets is 
suggested and given quantitative support. We argue that 
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FIG. 1: The fractional momenta of the partons at RHIC and 
LHC energies relevant for the production of dileptons through 
the NLO Compton and annihilation processes. The symbols 
on the curves denote the points when y — and yj ct =0.5 
(squares), 0.0 (crosses) and -0.5 (diamonds). 



these measurements could have advantages over measure- 
ments of real photons, even though they are similar in 
philosophy. In measurements of real photons, photons 
from the radiative decays of mesons, predominantly 7r°'s, 
constitute a huge background which overwhelms the real 
photon yield except at very high transverse momenta. 
The real photon sample can be enhanced by rejecting 
photons for which an accompanying photon is measured 
with the two photons having a photon-pair mass which 
falls within the 7r° mass window. However, this method 
fails at low transverse momentum for central heavy ion 
collisions due to the need for a large acceptance in order 
to measure both photons with high probability, and due 
to the high it multiplicity and the corresponding large 
number of combinatorial photon-pairs which results in 
the rejection of essentially all photons as possible de- 
cay photons. As the transverse momentum of the ir 
increases, such as when it is the leading hadron in the 
jet, the opening angle between the two decay photons 
decreases, making it possible to identify and reject the 
decay photons with high probability and with fewer false 
combinatorial possibilities. However, for photon mea- 
surements using calorimeters, at large transverse momen- 
tum the opening angle becomes so small that the sym- 
metric decay photons will merge into a single shower clus- 
ter and appear as a single isolated photon, making it im- 
possible to isolate them and discriminate between a single 
real photon and two overlapping photons from the decay 
of the 7T° . This then leaves a narrow window of photon 
(jet) energies where the real photons can be identified 
and such studies may be conducted. Of course, quantita- 
tive limits will depend on actual detector characteristics 
of geometry and granularity, and on the particle multi- 
plicity. Typical calorimeter geometries have 7r°-photon 
shower merging limits as low as 20 GeV/c. 



Now consider the case of a dilepton recoiling against 
the jet. The £ + and £~ are easily separated for arbitrarily 
large momenta. There are no other sources of dileptons 
having large qx (see later). Furthermore, the mass of 
the lepton-pair M gives an additional handle on the ini- 
tial state scattering [5j, as we can study dileptons which 
have the same transverse momentum (i.e. the recoiling 
jets have the same transverse momentum) but different 
masses. One pays the price of a low counting rate, but if 
the signal can be separated from the background and if 
luminosities are sufficient, we have a versatile tag. 

To repeat, the dileptons produced through the pro- 
cess qq — > £ + £~ will have only a modest qx resulting 
from the intrinsic transverse momentum of the partons. 
The correlated charm and bottom decays, which offer a 
huge contribution to the dilepton production, will also 
be governed by the same order of qx given from the in- 
trinsic momenta of the partons, in lowest order. Even 
though the cc and the bb will have a vanishing qx in the 
lowest order, the random transverse momenta of the lep- 
tons in the semi-leptonic decay of the resulting D and B 
mesons could still result in a net qr for the lepton-pair. 
At the NLO, the heavy-quark pairs will have additional 
qx, which would then get translated into the qr of the 
lepton-pair. One can see that the heavy quark will trans- 
fer only a part of its momentum to the lepton, (the semi- 
lcptonic decays B — > Y)£v and D — > K£v have three-body 
final states) and thus there is hope that for large qx of 
the pair, the contribution of the heavy-quark decays may 
become smaller than the NLO Drell-Yan contribution. It 
is also likely that one may be able to suppress or account 
for these contributions explicitly if one can identify the 
vertex of the decaying D or B meson, or reject leptons 
within a region of jet-like activity. Moreover, if the charm 
and the bottom quarks lose energy due to collisions or ra- 
diations in the plasma 0, their contribution will again 
be reduced [g, as the individual momenta of the quarks 
will be reduced, leading to a reduction in the net qx of 
the resulting lepton-pair. 

These hopeful considerations are quantitatively tested 
in the following. Our goal here is to perform simple phe- 
nomenological estimates, and to identify kinematical do- 
mains appropriate to RHIC and to the LHC where our 
tag would shine through the background. Note in passing 
that the CMS experiment has considered the possibility 
of observing jets tagged by W and Z bosons, which are 
governed by the same criteria M, and are thus open to 
similar vulnerabilities. See also Hoi. 



II. DILEPTON PRODUCTION 

At the outset, we add that the thermal radiation of 
dileptons from the QGP and the hot hadronic matter, 
a subject of considerable research [j| will not be im- 
portant [TH in the domain of very large qx of interest 
here. The dileptons from the annihilation of jets passing 
through the QGP [13 , could however have some contri- 
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bution, at least at the LHC. We shall report on this in 
a future publication. In the following we consider only 
the Drell-Yan process and the correlated decay of heavy 
mesons. 



A. Drell-Yan : Lowest Order 

The cross section for the lowest-order Drell-Yan pro- 
cess is given by 



da 



Ana 1 



dM 2 dy 9M 4 



F(X!,X 2 ,M 2 



(2) 



where 



F{x u x 2 ,M 2 



XiX 2 



E 



e\ [qf( Xl ,M 2 )qf(x 2l M 2 ) 

+ qf{x u M 2 )q?{x 2 ,M 2 )]{?>) 

qi are quark structure functions, and the sum runs over 
the quark flavours. We further have, 



x 2 = Me' 



(4) 



In the above s is the square of the centre of mass en- 
ergy per nucleon. The parton cross sections are scaled 
by the nuclear thickness function Tab(J> = 0) to obtain 
dN/dM 2 dy for the Drell-Yan production of dileptons in 
central (impact parameter b=0) A+B nuclear collisions. 
As observed in the introduction, if the partons do not 
have any intrinsic momentum then the net qT of these 
lepton pairs would be identically zero. The effect of in- 
cluding the intrinsic transverse momentum of the par- 
tons, to the LO Drell-Yan can be seen through the stan- 
dard practice of folding-in a Gaussian distribution in qT 
with a width of ~ 400 MeV. In the following we use the 
CTEQ4M parametrization [l| of the parton distribu- 
tions. 



B. NLO Drell-Yan : Annihilation 

The NLO process for dilepton production has several 
interesting features relevant to the present study which 
we recall here. We reiterate that we are interested only in 
large processes, and thus the treatment discussed here 
is sufficient for this purpose Q. A full treatment 
which regularizes the behaviour of the cross section as 
qx — ► is of course available. First consider the anni- 
hilation process, q + q — > g + 7*. The differential cross 
section for the production of a dilepton having mass M, 
transverse momentum qx, rapidity y, with a jet associ- 
ated at the rapidity yj ct through the annihilation process 
is given by 



dtr* 



dM 2 d 2 q T dy dy jct 



F(xi,x 2 ,fi F 



1 da J 



7T dM 2 dt 



(5) 



o 



Au+Au@Vs=200 AGeV 
M=2 GeV, q T =10 GeV 



y-0 




FIG. 2: Rapidity correlation of the dilepton and the recoiling 
jet at RHIC. 



We have Q 

dd A _ 4 2a 2 a s (n R ) (t - M 2 ) 2 + {u- M 2 ) 2 



dM 2 dt 
where 



9 ZM 2 



M 

M 2 = 



X\X 2 S 

1 



s 2 tu 



M 



(6) 



-- s x T xie 



u - M 
t u 



— — s XTX 2 e" 
s q% 



(7) 



and 



xi 



x 2 



!- v 
2 



= 2 XTe V 



vs<* 



2 XT€ 



2 XT€ 



-Vict 



(8) 



Above, any quantity under a caret (e.g. a) is asso- 
ciated with a parton-parton process. We set the renor- 
malization and the factorization scales to [1r = [If = 
M 2 + q T , as appropriate for the study of the transverse 
momentum of the lepton-pairs. It is obvious that if one 
decides to choose = = M 2 , the production of 
NLO Drell-Yan dileptons will change. We further have 



—2 



4t , xt — 



2q T 



M 2 



(9) 



We also add that the momenta of the lepton pair q, the 
incoming partons (pi), and the recoiling parton (A;), in 
the nucleon- nucleon centre of mass system are given by: 

Pi = zi-^(l,0,0,+l), 

P2 = X2~Vs(l, 0,0,-1), 

q = (?o,qT)9i): 

k = (fco,-qr,fei)- (10) 
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FIG. 3: Rapidity correlation of the dilepton and the recoiling 
jet at the LHC. 



FIG. 5: Jet-rapidity integrated dilepton rates at the LHC for 
various masses of the pair. 
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FIG. 4: Jet-rapidity integrated dilepton rates at RHIC for 
various masses of the pair. 



C. NLO Drell-Yan : Compton 

The expression for the Compton contribution to the 
Drell-Yan process at NLO is a little more involved. We 
have 3 

da c _ . ,1 dd c (s,u) 

= F 1 (xi 1 x 2 )-- 



dM 2 d 2 qr dy dyj 



F 2 (xi 1 x 2 ) 



it dM 2 du 

1 dd c {s,i) 
ir dM 2 du ' 



(11) 



where 



F 1 (x 1 ,x 2 ) = x 1 x 2 g A {x 1 )^2e 2 [qf (x 2 ) + qf (x 2 )] , 

i 

F 2 {x u x 2 ) = X!X 2 ^2e 2 [qf(xi) +qf(xi)] g B {x 2 ) , 

i 

(12) 



dd c (s,u) _ 1 2a 2 a s (s - M 2 ) 2 + (u- M 2 ) 2 



dM 2 du 6 3M 2 



s 3 (— u) 



; (13) 



and 



da c (s, t) _ 1 2a 2 a s (s - M 2 ) 2 + (t - M 2 ) 



dM 2 du 6 3M 2 



S3(- t ) 



(14) 



In the above, we have not explicitly shown the factoriza- 
tion and the renormalization scales, which are taken as 
Pf = M-R = + It as l° r l ne NLO annihilation term. 



D. Heavy-Quark Production 

We consider the production of heavy quarks at NLO. 
The lowest order contributions come from the processes 
qq — > QQ, gg — > QQ. The NLO terms originate from 
qq -> QQg, gq(q) -> QQq(q), and gg -> QQg. The 
general expression for the production can be written as 

d<y = Y.J dx i dx 2 dd l0 {s,M 2 Q ,Q 2 )f?( Xl ,Q 2 )fZ(x 2 ,Q 2 ) 

(15) 

where are the distribution functions for the partons 
in the nucleon in the nucleus A, the fractional 

momenta of the incoming partons, and s = x\x 2 s is the 
parton-parton centre of mass energy. The parton cross- 
section &i ; j(s, Mq, Q 2 ) can be written as: 



2 ^_«2(Q 2 ), ( Q 2 

Ji,j Pi 



where 



4M% 

X\X 2 S 



(16) 



(17) 
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and 



f itj {p,Q 2 /M*) = /ffOO + W? 2 ) 4»+ 



/g(p)ln(Q a /M5) 



(18) 



are taken from Ref. |l6j|. We have used the factorization 
and the renormalization scales as Q 2 — Mq+ < p T >, 
where < p T > is the average of the transverse momen- 
tum of the produced quark and anti-quark. The corre- 
lated semi-leptonic decay of the heavy (D or B) mesons 
is then estimated by fragmentation of the heavy-quarks 
to produce the meson and then estimating its decay in 
the quark model ^t|- The branching ratios to leptons 
for the decay of the D and B mesons are taken as 12% 
and 10% respectively. Even though more sophisti- 

cated approaches exist that outlined above is suffi- 
cient for our purpose, as we only wish to know where 
these contributions are small. We have verified that 
our results are quantitatively similar to those obtained 
from PYTHIA 19] and also by Gavin et al. [H|, where 
parametrizations of the measured lepton spectra from 
semi-leptonic decay of B and D mesons are used to gen- 
erate the lepton momenta. The total cross sections for 
heavy qq we thus obtain are shown in Table |U These val- 
ues do depend on the choice of structure functions and 
scales 201. 



(GeV) 




NLO 


^bb 


NLO 
^bb 


200 


59.3 


145.9 


0.84 


1.67 


5500 


1054 


3362 


79.7 


174.4 



TABLE I: Total cross sections for heavy qq production in pp 
collisions. All values are in /ib. 



> 

O 



10° 






Au+Au@ 200 AGeV 
b=0 fm 


10" 2 








lO- 4 








10 -6 




LO DY x 1.5^ 





10 

M (GeV) 



15 



20 
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tons from their correlated decay at RHIC energies. The LO 
Drell-Yan contribution is also shown. 
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Finally, the momentum lost by the heavy quarks dur- 
ing fragmentation is taken to be negligible |2l| . 



FIG. 7: Same as FigElfor the LHC 



III. RESULTS 

As a first step we explore the range of fractional mo- 
menta of the partons (Eq. |HJ which are relevant for our 
purpose (see Fig. 1). The diagonal provides the values 
for the situation when the rapidity of the dilepton (y) 
as well as that of the recoiling jet (i/jet) is zero. The 
squares denote the point when the yj et = +0.5, while the 
point y^t = —0.5 is indicated by the diamonds, when the 
dilepton y = 0. Thus large M and very large qr of the 
lepton pair do not sample very small values of x, and the 
nuclear modification of the parton distribution functions 
will hence not play a large role. 

Next we study the rapidity correlation between the 
dilepton and the jet. We see (Figs. 2 and 3) that there is 
a strong positive correlation between the rapidity of the 
dilepton and the recoiling jet; i.e. when the dilepton has 
a positive rapidity the jet also has a positive rapidity, and 



the correlation has a width of about 2 units of rapidity. 
This aspect has been discussed in detail by Kajantie et 
al. and has been attributed to the dominance of the 
Compton term through the 1 /s term in the scattering 
cross section. Note that the finite y values are chosen as 
the mid-points of the forward/backward rapidity cover- 
age of PHENIX at RHIC, and they fall within those for 
the ALICE and CMS experiments at the LHC. 

We now study the jet-rapidity integrated dilepton rates 
at RHIC (Fig. 4) and at the LHC (Fig. 5) due to the NLO 
processes discussed above. We have also given the results 
for the LO Drell-Yan contribution, supplemented with a 
modest intrinsic kx for the partons. As expected, the 
NLO processes dominate as the transverse momentum of 
the pair (which is also the transverse momentum of the 
recoiling jet) increases. The NLO contributions to the 
Drell-Yan and also the heavy-quark production (see be- 
low) have been calculated ignoring the intrinsic momenta 
of the partons. Including it will increase the < qt > of 



6 



,,,,,,,, 


, | , , , , | , , 

Au+Au® 200 AGeV 






1.5 < M < 2.5 GeV 






Central Collisions 


- 


\j-ps ->i + r+x 

- 




- 




DY@NL0 




bb ->i + r+x^Ap-, 








t-^ 1 n i_ i i i i i 





20 

q T (GeV) 



30 



I ' ' ' ' I ' ' ' 

Au+Au® 200 AGeV 
3.5 < M < 4.5 GeV 
Central Collisions 




q T (GeV) 

FIG. 8: Rate of dilepton production with pair mass 1.5 < 
M < 2.5 GeV (upper panel) and 3.5 < M < 4.5 GeV (lower 
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energies. 



the resulting lepton-pair and may enhance the yield by 
50-100% depending on the value of the intrinsic momen- 
tum. 

The semi-leptonic decay of D and B mesons gives a 
large contribution to the dilepton production We 
show the results of our calculations, which extend up to 
large lepton pair masses. The LO Drell-Yan contribution 
has been scaled by a factor of 1.5 to account for the NLO 
terms as our exploratory calculations will be inadequate 
at low qx, and a momentum integral has to be performed 
in order to obtain the mass distributions. 

Figs. El and [7| suggest that it may be very difficult 
to observe the dileptons recoiling against a jet in the 
background of the correlated decay of charm and bot- 
tom. This indeed is the case for dileptons radiated from 
thermal sources |l l|. However, this conclusion changes 
when the qr distribution of the dileptons is studied for 
large pair masses, as shown in Fig.[S]for mass windows of 
2 and 4 GeV at RHIC energies. We find that dileptons 
originating from the NLO Drell-Yan process dominate at 
large qx, and with anticipated total sampled event sizes 
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FIG. 9: Same as Figilfor the LHC. 



of ~ 10 10 events, open a window for a clean measurement 
of dilepton-tagged-jets. This window shifts to higher qx 
at the LHC (Fig©, and identificat ion becomes more dif- 
ficult than at RHIC. The situation may improve if the 
heavy-quarks lose energy in the medium, degrading the 
D and B meson decay contributions to lower dilepton qx ■ 
Furthermore, the background can be reduced if the con- 
tributions of the heavy- meson decays can be isolated and 
rejected, such as by determination of the decay vertices, 
or by jet-like isolation cuts. 



IV. CONCLUSIONS 

In brief, we have estimated the production of jets re- 
coiling against a dilepton produced in the Drell-Yan pro- 
cess at NLO in relativistic heavy ion collisions at RHIC 
and LHC energies. The dileptons could be used to tag 
the jet and give the energy of the progenitor parton, 
which would provide precise information with which to 
determine the rate of energy loss of partons during their 
passage through quark gluon plasma. The background 
from the correlated decay of charm and bottom mesons, 
which plagues the identification of thermal dileptons > 
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is found to be unimportant for large jet (or dilepton) 
transverse momenta, that is, in the kinematical region 
of interest for jet-quenching studies. Even though the 
counting rates will be small, the anticipated event sam- 
ples suggest that dilepton jet-tagging is feasible. This 
method could thus prove useful in the resolution of sev- 
eral issues of fundamental importance in the physics of 
high energy partons traversing a plasma of quarks and 
gluons. 
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